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Abstract 


One of the dangers facing Iran’s plains is subsidence. Land subsidence due to various reasons, 
such as over-harvesting of underground water sources and climate changes, has caused many 
problems in agricultural lands, roads, and power and energy transmission lines. Therefore, it is 
important to deal with the causes and influencing factors to control the risks. Ardabil plain is also 
one of the areas prone to subsidence due to climate changes and the decrease of underground 
water in recent years. The aim of the present research is the zoning subsidence risk in this plain. 
First, the factors influencing subsidence in the Ardabil plain (slope, land use, lithology, distance 
from the fault, distance from the waterway, drop in water level, distance from the city and the 
village) were identified. Then the layers of data were applied in the geographic information 
system. In the next stage, the weighting of the investigated factors was done using ANP method 
and in the Super Decision software, and the final analysis and modeling were done using the 
MABAC method as one of the multi-criteria decision making methods. Finally, the resulting map 
was Classified into five categories from very low risk to very high risk. The results showed that 
such factors as water level drop, distance from the river and lithology have the highest weighting 
factor. Moreover, 244.29 and 370.59 square kilometers, respectively, of the area of this plain, are 
in very dangerous and dangerous classes. Finally, due to the high potential of Ardabil plain in 
terms of subsidence, protective and management measures should be taken into consideration by 
relevant authorities and institutions. 
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جغرافيا و مخاطرات محيطى. سال يازدهم. شمارة چهل و چهارم زمستان ۶۰۱ صص PW‏ 
مقاله پژوهشی 


ویژه‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران يا بحران مدیریت) 


ارزيابى و بهنه‌بندی خطر فرونشست با استفاده از الكوريتم تطبيقى MABAC‏ و ANP‏ 
(مطالعة موردى: دشت اردبيل) 


D‏ موسی عابدينى '- استاد ژئومورفولوژی» دانشكده علوم اجتماعی» دانشكاه محقق اردبيلى؛ اردبیل» ايران. 
ليلا آقايارى- دانشجوى دكترى ژئومورفولوژی. دانشكده علوم اجتماعی, دانشگاه محقق اردبیلی» اردبیل» ايران. 
صیاد اصغری سراسکانرود- استاد ژئومورفولوژی دانشکده علوم اجتماعی, دانشگاه محقق اردبیلی» اردبیل ايران. 


تاريخ دریافت: ۰۰/۱۰/۲۶ع ۱‏ تاريخ بازنگری: ۱۶۰۱/۳/۸۱ تاريخ تصویب: ۱۶۰۱/۶/۶ 


حکیده 


یکی از مخاطرات پیشروی دشت‌های کشور مخاطره فرونشست است. پدیده فرونشست زمين به دلایل 
مختلف ازجمله برداشت بیش ازحد منابع آب زیرزمینی و تغییرات جوی سبب بروز مشکلات و معضلات 
فراوان در زمین‌های كشاورزىء جاده‌هاء خطوط انتقال نیرو و انرژی شده است؛ ازاین‌رو پرداختن به علل و 
عوامل SE‏ گذار جهت کنترل و مدیریت خطر دارای اهمیت است. دشت اردبیل نيز به خاطر تغییرات 
اقلیمی و افت سطح آب‌های زیرزمینی در سال‌های اخیر یکی از مناطق مستعد جهت مخاطره فرونشست 
است. هدف تحقیق حاضرء پهنه‌بندی خطر فرونشست در اين دشت است. در اين مطالعه dal‏ عوامل موثر 
جهت ایجاد فرونشست در دشت اردبیل (شیب. کاربری اراضی, لیتولوژی. فاصله از گسل. فاصله از 


آبراهه. افت سطح آب. فاصله‌ازشهر و روستا)» شناسایی شدند و سپس نسبت به Og‏ لایه‌های اطلاعاتی در 


سامانه اطلاعات جغرافيايى اقدام كرديد. در مرحله بعد وزن‌دهی عوامل موردبررسى. با استفاده از روش 
ANP‏ و در محبط نرم‌افزار Super Decision‏ انجام گردید و تحليل و مدلسازى نهايى با استفاده از 


Email: abedini@uma.ac.ir ۰۹۱۳۵۵۵۲۲7 نويسنده مسئول:‎ ١ 
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روش MABAC‏ به‌عنوان یکی از روش‌های تصمیم‌گیری چندمعیاره. انجام شد. درنهایت. نقشۀ حاصله 
در ينج رده با خطر بسیار کم تا حطر بسیار زياد طبقه‌بندی گردید. با توجه به نتایج مطالعه عوامل افت 
سطح آب. فاصله از رودخانه و لیتولوژی بیشترین ضریب وزنی را به خود احتصاص دادند. همچنین, quU‏ 
مطالعه OUS‏ داد؛ به ترتیب ۲۶/۲۹ و ۹ کیلومترمربع از مساحت اين دشت. در طبقات بسیار پرخطر 
و پرخطر قرار دارند. درنهایت می‌توان اظهار داشت. نظر به توان بالای دشت اردبیل» از لحاظ رخداد 
فرونشست. بایستی اقدامات حفاظتی» مدیریتی در سطح دشت اردبیل مورد توجه مسئولان و دستگاه‌های 
ذیربط قرار گیرد. 

کلیدواژه‌ها: مخاطرات. پدیده فرونشست. MCDM‏ دشت اردبیل. 


۱- مقدمه 

چشم‌انداز سطح زمین به‌عنوان سیستمی پیچیده حاصل تعامل فاکتورهایی نظیر فرآیندهای جغرافیایی» اقلیمی؛ 
زمانی و فعالیت‌های انسانی است. رشد کالبدی, بدون رعايت اصول برنامه‌ریزی شهری و ساخت‌وسازهایی بدون در 
نظر گرفتن پتانسیل مخاطرات موجود در پهنه‌های در معرض مخاطرات quee‏ آسیب‌پذیری و خسارت‌ها در 
محدوده‌های مورد نظر را تشدید می‌کند. عواملی که با تغيير در چگونگی ساختار و کارکرد توي وكرافى طبیعی و ريز 
اقلیم به‌ویژه در محيطهاى انسان‌ساخت. زمینه‌ساز ایجاد مخاطراتی نظیر زمین‌لغزش» سیلاب‌شهری» زمین‌لرزه و 
فرونشست می گردد (فنی و همکاران. ۱۳۹۱). براساس تعریف سازمان زمین‌شناسی ایالت متحده پدیده فرونشست 
زمين شامل فروریزش يا نشست روبه پائین سطح زمین است که می‌تواند دارای حرکت قائم رو به پائین سطح زمین 
و کمی افقی باشد (عابدینی»۱۳۹۱). پدیده فرونشست زمين به دلایل مختلف ازجمله برداشت بیش از حد منابع آب 


زیرزمینی و تغییرات جوی سبب بروز مشکلات و معضلات فراوان در زمین‌های کشاورزی, جاده‌ها؛ خطوط انتقال 


سرعت شکل‌گیری» به دو دسته سریع و تدریجی تقسیم می‌گردد. حالت تدریجی به دلیل تلفات انسانی کم آن در 
مقایسه با دیگر پدیده‌های طبیعی مانند ;3« لغزش و سيلاب كمتر موردتوجه قرار فى كيرد (مدنی و همکاران. 
٦‏ ( ۱۵ 


اما فرونشست در حالت سریع باعث آسیب به سازه‌ها و تلفات انسانی می‌گردد. فرونشست زمين در سال‌های 


اخير به مشکلی جهانی تبدیل شده است. چراکه در اكثر کشورها در حال وقوع است. فرونشست سبب تخریب 
کانال‌های آبیاری» شکست خطوط لوله‌های col‏ و نفت. خطوط انتقال نیرو و تخریب سازه‌ها و تشدید فرایند زمين 
لغزه و سيلابها مى شود (محمودى و همکاران. KITAY‏ فرايندهاى طبيعى مانند حركات تكتونيكى» فوران dx sS‏ 


سال يازدهم ارزيابى و يهنهبندى خطر فرونشست با استفاده از .... £0 


وسيع م ىكردند (ژویی ' و همکاران. (THT‏ به‌طورکلی فرونشست متأثر از فعالیت‌های انسانى مانند برداشت بىرويه 
آب از سفره‌های آب زيرزمينى و پارامترهای زمین‌شناسی است كه در صورت عدم شناخت اين يديده سبب آسیب 
جدى به زيرساختهاى اقتصادى م ىكردد (گالووی" و همكاران. ۲۰۰۹). پژوهشگران رابطه بين توسعه شهرى و 


tij‏ زمين را رابطه‌ای دو جانبه مىدانند كه هردو بر هم تأثير كذاشته و تأثير مىيذيرند. به كفته آن‌ها تمدن 


شهرها به جهت رشد سريع در توسعه شهرى در بخش‌های مختلف صنعت. تجارت و حمل ونقل اثرات 
زیست‌محیطی منفى از خود بر جا می‌گذارند؛ مانند: تغيير در كاربرى اراضى كشاورزى به مسکونی» مصرف بىرويه 
منابع آب‌های زيرزمينى جهت كسترش فعالیت‌های صنعتی و افزايش جمعيت اين اثرات منفى درنهايت منجر به بروز 
يديده فرونشست می‌شود و همین يديده می‌تواند برنامه و فرایند توسعه شهرى را متأثر سازد (منتظریون و همکاران؛ 
۸ 

به‌منظور شناسایی و کاهش پیامدهای ناشی از پدیده فرونشست نیاز به یک سیستم يايش فرونشست بیش‌ازپیش 
در کشور و دشت اردبیل احساس می‌شود. اولین گام در پایش فرونشست اندازه‌گیری مقدار جابه‌جایی ناشی از آن در 
سطح زمین است. با علم به میزان نرخ فرونشست و همچنین پهنه و گستره مناطق تحت تأثير, نه‌تنها قادر خواهیم بود 
میزان پیشرفت اين پدیده را در یک منطقه شناسایی کنیم و از پیشرفت بیشتر آن جلوگیری نماييم. بلکه می‌توانیم 
اطلاعات باارزشی از ویژگی‌های زمین‌شناسی و سفره آب زیرزمینی به دست آوریم. 

امروزه؛ از یک سوء سیستم اطلاعات جغرافیایی. ابزاری مؤثر برای سیستم‌های پشتیبانی تصمیم‌گیری تعاملی برای 
عملیات مدیریت ريسك هستند (بوخیر " و همکاران» (TeV‏ و در سال‌های اخير. جهت مطالعه فرونشست مورد 
استفاده فراوان قرار می‌گیرند. از سوی دیگر برای كشف طیف وسیعی از گزینه‌ها از نظر درگیری‌های عینی و 
معیارهای چندگانه. از روش MCDM‏ استفاده می‌شود (آهر * و همکاران» ۲۰۱۳). در اين روش‌های تصميم كيرى 
برای بهینهسازی مدلء از معیارهای مختلفی برای افزايش صحت تصمیمات استفاده می‌شود (جورجیو و همکاران" 
۵ عرب عامری و همکاران: ۲۰۱۸). لذا می‌توان گفت» استفاده از سیستم‌های اطلاعات جغرافیایی و روش‌های 
تصمیم گیری چندمعیاره با یک رویکرد تلفیقی, می‌تواند باعث تسریع در روند برنامه‌ریزی در تشخیص موارد بحرانی 
و اضطراری شده و منجر به صدور نتایج مناسبی گردد. مطالعات متعددی در نقاط مختلف جهان انجام شده و توانایی 
تحلیل‌های تصمیم‌گیری چندمعیاره» برای ارزیابی و پهنه‌بندی خطر فرونشست نشان داده است. تاکنون مطالعات 


فراوانی پیرامون پهنه‌بندی مخاطرات ازجمله فرونشست در ايران و ساير کشورها انجام شده است که به برخى از laol‏ 


1 zhu et all 

2 Galloway at all 
3 Bukheir et all 
4 Aher et all 

5 Georgiou et all 


£a‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


اشاره می‌شود. استیفن ۲ QVO‏ به بررسى علل فرونشست و عوامل تشديدكننده و خسارات ناشی از فرونشست را 
بررسى كرده و مهمترين علت فرونشست را تخليه پیش‌ازحد آب‌های زيرزمينى و انحلال سنگ‌های کربناته اعلام كرد. 
كوبى انجوين (V‏ در كتاب fle‏ فرونشست زمين در شهر هوشىمين به مطالعه اين يديده يرداخته و توصیه‌های 
مناسبى در جهت مديريت آن در سطح شهر بیان كرده است. گامبولاتی و همكاران (۲۰۰7) به مطالعه تأثير 
فعالیت‌های انسانی بر ميزان فرونش ست زمين در مناطق مسكونى پرداختند و طبق بررسىهاى انجام شده برداشت 


بىرويه منابع زيرزمينى و تجمع سازه‌ها را مهمترين علت فرونة - در مناطق مسكونى بیان كردند. هوانگ و 


همکاران " (۲۰۱) به بررسى فرونشست شهر دزو با روش SBAS‏ پرداختند. بررسىها طى بازه زمانى دو ساله نشان 


می‌دهد كه مركز فرو بخش‌های شرقى و شمالشرقى با ميزان فرونشست £0 ميلىمترى است و اين نشست 
زمين و افت سطح آب‌های زيرزمينى در شهرستان كرج با استفاده از روش تلفيق وزنی در محيط 5 كرده و نتيجه 
كارشان نشان داده كه بين وضعيت توپوگرافی. ضخامت سازند و برداشت آب ارتباط تنگاتنگی وجود دارد و بيشترين 


ميزان فرو در مناطق با برداشت obj‏ آبرفت‌های ضخیم» و مناطق يست دشت مشاهده می‌شود. عابدینی و 
همکاران (۱۳۹۵) در مورد فرونشست دشت روانسرا- سنجابی با استفاده از مدل تحلیل سلسله مراتبی کار کرده‌اند و 
نتایج نشان داده که ارتباط و همبستگی قوی ميان عمق رسوبات ریز دانه چاه‌ها و احتمال فرونشست زياد آتی وجود 
ندارد به‌طوری که حداکثر عمق رسوبات ريز دانه در نواحی مستعد فرونشینی با احتمال خطر UL‏ حدود ٩۱۲‏ متر است 
اين در حالى است که در برخی از نقاط آبخوان که استعداد کمتری برای فرونشست دارند عمق رسوبات ريز دانه بسیار 
RAS‏ و حدود ۲۰ ۲۲ متر است. 

شادفر و همکاران (TAL)‏ در مقاله‌ای به پهنه‌بندی خطر فرونشست زمين با استفاده از روش AHP‏ يرداختند. 
نتایج حاصل نشان می‌دهد که فاکتور برداشت بی‌رویه با ضریب 4۸ درصد در درجه اول و عامل سنگ‌شناسی با 
ضریب VE‏ درصد در درجه دوم اهميت در ایجاد فرونشست منطقه است. همچنین نتایج نشان داد که ۵ درصد از 
مساحت منطقه در پهنه‌بندی با خطر نسبتاً زياد و خيلى زياد قرار دارند. شايان (۱۳۹۵) به پهنه‌بندی فرونشست زمین 
در حوضه آبخیز قره‌چای همدان پرداختند. جهت پهنه‌بندی از ۵ ايه ارتفاع» زمین‌شناسیء هیدرولوژیکی شبکه‌های 
زهکشی و شیب استفاده کردند. نتایج OLY‏ می‌دهد که تغییرات ناهمسان در داده‌های مورفومتری شبکه زهکشی 
مهم‌ترین فاکتور موثر در تشدید افت آب‌های زیرزمینی و وقوع فرونشست در منطقه بوده و عوامل لیتولوژیکی و 
انسانی به‌تنهایی نقشی ندارند. فنی و همکاران YAT)‏ به پهنه‌بندی و تحلیل فضایی حساسیت مخاطرات چهار گانه 


1 Stephen 

2 coein enjoein 
3 Gambolati 

4 hoang at all 


سال یازدهم ارزیابی و پهنه‌بندی daz‏ فرونشست با استفاده از .... ۷ 


رودخانه. کاربری اراضی و زمین‌شناسی استفاده کرده‌اند و نقشه حساسيت v‏ در قالب ۵ كلاس حساسیت 


خیلی‌زیاد» زياد daos ga‏ کم و خیلی کم طبقه‌بندی و به‌صورت نقشه و جدول استخراج و ارائه شده‌اند. عالی يور و 
همکاران OY)‏ به پهنه‌بندی ريسك فرونشست زمین در اثر افت سطح أب زیرزمینی با استفاده از مدل تحلیل 
سلسله مراتبی فازی دشت اردبیل پرداختند و اقدام به تهیه نقشه احتمال فرونشست زمین و تعیین مناطق دارای 
ريسك پرداختند. لایه‌های اطلاعاتى شامل ۵ لایه زمین‌شناسی» فاصله از كسلء افت سطحآب زیرزمینی» شيب و 
ضخامت آبرفت به‌عنوان معيارهاى مؤثر مورد استفاده قرار گرفتند و نتايج نشان داد كه ۱۳ درصد منطقه داراى ريسك 
خيلىبالا است و حدود W/0‏ درصد مناطق جنوبى در بخش ريسكهاى خیلی بالاى فرونشست. به دليل استفاده 
بي شازحد از آب‌های زيرزمينى است. بهيارى و همكاران (YS‏ به ارزيابى تأثير ساختارهاى فعال بر خطر 
فرونشست زمين با استفاده از مدلهاى تصميمكيرى جندمتغيره پرداختند. در اين مطالعه از روش AHP‏ به‌منظور 
تجزیه‌وتحلیل خطر فرونشست در منطقه موردمطالعه با استفاده از لایه‌های مطالعاتی ارتفاع, شیب. آبراهه‌ها و 


زمین‌شناسی, نقشه پهنه‌بندی das‏ فرونشست در منطقه تهیه گردید. اين نقشه OUS‏ می‌دهد كه خطر فرونشست در 
بخش جنوبی منطقه نسبت به ساير مناطق بیشتر است. كوه بنانی و همکاران (۱۳۹۸) به پهنه‌بندی گستره خطر 
فرونشست زمین با بهره‌گیری از تداحل‌سنجی راداری در دشت کاشمر و خلیل‌آباد با استفاده از تصاویر سنتینل ۱ 
موبوظ a‏ مال هاف aul eb AY UAE‏ نشان داف این dice‏ جر پیش رين pagi ca ۱۷ Gib‏ در سنال 


فرونشست دارد و پیشتر در حومه جنوب غرب شهر کاشمر و مناطق غربی دشت خليل أباد متمرکز است. منتظریون و 


دشت آبخوان نورآباد با استفاده از روش تداخل‌سنجی 
راداری پرداختند و جهت بررسی وضعیت افت سطح آب‌های زیرزمینی از روش درونیابی IDW‏ استفاده کرده‌اند و 
نقشه‌های حاصل از فرونشست نشان می‌دهد که نرخ نشست در هرسال به‌طور میانگین ‏ سانتی‌متر و طی بازه زمانی 
٤‏ ساله ۱١‏ سانتی‌متر است. ابراهیمی و همکاران (۱۳۹۹» به يايش میزان فرونشست محدوده شهری دشت پاکدشت 
با روش تداخل‌سنجی پرداختند. نتایج حاصل از روش تداحل‌سنجی راداری محدوده مطالعاتی در طی بازه زمانی 
یک‌ساله بين ۱۵ تا ۸۵ میلی‌متر فرونشست را OUS‏ می‌دهد و مهم‌ترین و موثرترین عامل در فرونشست را هم افت 
سطح آب‌های زیرزمینی عنوان کرده‌اند. شاه کرمی و همکاران (۱۳۹۹) مناطق فرونشست در اثر اضافه برداشت آب 
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زيرزمينى در دشت اراک را موردبررسى قرار دادند. نتايج به دست آمده Sly‏ اين موضوع هستند كه نواحى مرزى و 
حاشيهاى دشت اراك در معرض افت بیش‌تر آب قرار دارند و تمركز رسوبات ريزدانه رسى در قسمت‌های غربى 
به‌مراتب بيشتر از نواحى ديكر است. درنتيجه احتمال می‌رود اين نواحى در معرض خطر فرونشست بيشترى باشند. 


زارعى و همكاران (۱۳۹۹) به تعيين رابطه فرو 


زمين و افت سطح آب زيرزمينى با دو روش تداخل‌سنجی 
رادارى و ايستكاه سايت gps‏ دشت سلماس پرداختند و نتايج نشان داد كه به ازاى یک متر افت سطح آب زیرزمینی؛ 
۸ متر سطح زمين در دشت نشست می‌کند. جوادى Or)‏ در رساله خود به برآورد و تحليل ميزان فرونشست 
دشت همدان- بهار با استفاده از تصاوير راداری و سرى زمانى SBAS‏ در طى بازه زمانى ۵ ساله و ميزان فرونشست 
اين دشت را ۵۱ ميلىمتر به‌صورت سالانه برآورد كرده است و بيشترين ميزان فرونشست مربوط به مناطق داراى 


افت منابع آب زيرزمينى است. فرزينكيا و همكاران ( 28٠١‏ به يهنهبندى خطر فرونشست زمين در دشت جوين با 


استفاده از مدل تحليل شبکه‌ای فازى يرداختند در اين تحقيق از روش ANP‏ و فازی. فرونة ت بوصورت 
کلاس‌بندی براى نشان دادن ویژگی‌های مناطق بحرانی فرونشست و اولویت‌بندی مناطق مختلف از نظر استعداد برای 
وقوع پدیده فرونشست به دست آمد و نتایج بیانگر اين است که OV‏ درصد قسمت مرکزی دشت جوین با خطر زياد 
تا خیلی زياد بهحساب می‌آیند. 


هرچند که تحقیقات معدودی در مورد دشت اردبیلی انجام شده ولی موضوع آن‌ها متفاوت بوده اسست. 


YA) e‏ در مورد ple‏ فرونة 2 دشت اردبیل کار 02,5 و مهم‌ترین fle‏ أن را برداشت بی‌رویه آب‌های 
زیرزمینی» فعالیت‌های تکتونیک حاصل از گسل‌های حواشی دشت و سنگینی ناشی از انباشت آبرفت‌ها عنوان کرده 
است 
سسا 


با توجه به منابع بررسی شده محققان مختلف اين تکنیک را بهعنوان ابزار نيرومندى برای مطالعه پدیده‌های 
مختلف زمین‌شناسی همچون فرونشست مطرح كرده است. روش MABAC‏ كه براساس رويكرد برنامه‌ریزی 
مطلوب و بهینه‌سازی چندمعیاره پایه‌ریزی شده است. به‌عنوان روش تصمیم‌گیری چند معیاره در تلفیق با فرایند 
تحلیل شبکه به‌منظور پهنه‌بندی پتانسیل وقوع فرونشست دشت اردبیل استفاده شد. طبق بررسی als‏ منتشر شده 


تاکنون در دشت اردبیل با اين روش و داده‌های ورودی» پژوهشی ارائه نشده است لذا جهت برآورد دقیق‌تر میزان 


فرو و عوامل تأثیرگذار بر آن لزوم و خلأ انجام تحقیق با این روش كاملاً احساس می‌شود. لذا اهداف تحقیق 
حاضر بررسی فرونشست دشت اردبیل با استفاده از تکنیک تداخل سنجی راداری, پهنه‌بندی فرونشست. اندازه‌گیری 


معیارهای مؤثر بر فرونشست و بررسی مژثرترین عوامل در فرونشست دشت اردبیل است. 
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۲- مواد و روش‌ها 
۱-۲- منطقه موردمطالعه 


اردبیل واقع شده است. در بررسی زمین‌شناسی دشت اردبیل با مختصات جغرافیایی EA‏ درجه تا ÉA‏ درجه و ۳۰ 


دقيقه طول خاوری و ۳۸ درجه و ۲۰ دقيقه شمالی قرار دارد. زمین‌های حاشیه دشت از اطراف به سمت دشت شیب 
دارند. ارتفاع متوسط اراضی کشاورزی حدود ۱۳۲۰ متر از سطح دریاهای آزاد قرار دارد. وسعت آن حدود ۸۲۰ 
کیلومترمربع است و جزئی از حوضه آبخیز رودخانه قره‌سو به شمار می‌رود. دشت اردبیل متشکل از ۸۸ واحد 
روستایی ده واحد دهستان و ۲ واحد شهری اردبیل و آبی بیگلو است. شهر اردبیل مرکز استان اردبیل در غرب اين 
دشت ally‏ شده است. در محدوده دشت اردبیل تعداد ۲۲۵۳ حلقه چاه عمیق و نیمه‌عمیق, ۲۰ دهنه چشمه و VA‏ 


رشته قنات وجود دارد كه به مصارف کشاورزی اختصاص دارند (رستمزاده و ATE iS‏ 
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شكل -١‏ نقشه موقعيت دشت اردبيل 
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-Y-Y‏ داده‌ها و ابزارهاى مورد استفاده 

در راستای انجام اين پژوهش از داده و تصاویر ماهواره‌ای و نرم‌افزارهای متعددی استفاده شده است که مشتمل بر 
موارد زیر خواهد بود: 

پژوهش حاضر از نوع کاربردی بوده و روش تحقیق آن تحلیلی مبتنی بر آنالیز داده‌ه؛ سیستم اطلاعات جغرافیایی و 
استفاده از فنون تحلیل چندمعیاره اسست. جهت پردازش تصاویر و تجزیه‌وتحلیل داده‌ها نیز از نرم‌افزارهای 
Excel ; Arc GIS Ecognition ENVI‏ بهره‌گیری شده است. جهت بررسی خطر فرونشست. در مرحله جمع‌آوری 
دادههاء ابتدا بايد ملفه‌های تأثیرگذار مشخص شده و بر مبنای آن‌ها داده‌ها و اطلاعات مورد نیاز گردآوری و طبقه‌بندی 
شوند. بدین منظورء در اين مطالعه؛ ابتدا عوامل مؤثر (شامل: شیب لیتولوژی, کاربری اراضی فاصله از آبراهه» فاصله 
از كسلء افت سطحآب و فاصله از شهر و روستا) با توجه به شرایط طبیعی و انسانی منطقه شناسایی گردید. در 
مرحله بعد لایه‌های اطلاعاتی مربوط به هر یک از عوامل, در محیط سيستم اطلاعات جغرافیایی تهیه گردید. LU‏ 
اطلاعاتی شبکه آبراهه با استفاده از نقشة رودخانه‌های استان اردبیل» تهیه گردید. لایه‌های شیب نیز با استفاده از مدل 
رقومی ارتفاعی, دانلود شده از سايت (USGS)‏ تهیه شد. لایه‌های اطلاعاتی مربوط به لیتولوژی (مقاومت سنگ‌ها) و 
كسلهاء با رقومی‌سازی از روی نقشه زمین‌شناسی اردبیل؛ با مقیاس ۱:۱۰۰۰۰۰ تهیه شد. برای استخراج کاربری 
اراضی محدوده موردمطالعه» از تصاویر سنجنده لندست هشت (سنجنده (OLI‏ به شماره گذر VW‏ و ردیف ۳۶ 
تاريخ ۰ ۷ dl Y‏ شده از «(earthexplorer.usgs.gov ) cole‏ استفاده شده است (حدول ۱). 

تصاوير ماهوارهاى ايتيكال: تصوير ماهوارهاى ايتيكال در اين تحقيق بهمنظور استخراج نقشه كاربرى اراضى سال 
۰ استفاده شد. مشخصات تصاوير لندست اخذ شده از سايت <(earthexplorer.usgs.gov)‏ در جدول Y‏ ارائه 
شده است. به‌منظور تصحيحات اتمسفرى و راديومتريكى از نرمافزار ENVIS.3‏ و به‌منظور استخراج نقشه‌های مربوطه 
از نرم‌افزار ARCGIS10.5‏ استفاده شد. برای طبقهبندى کاربری اراضى از روش طبقهبندى شىعكرا در 


نرم‌افزار developer64‏ 0 استفاده شد. همجنين بهمنظور به دست آوردن صحت طبقه‌بندی از Google‏ 


Earth‏ استفاده شد. 
حدول ۱- اطلاعات تصاویر سنحنده لندست ۸ 
SCENE-CENTER-TIME SUN-AZIMUTH ROW PATH Sansor name Year‏ 
VT \£0,0+ ۷ v vw Landsat 8 ۳۹۰‏ 


در روش طبقه‌بندی شیءگرا؛ اطلاعات طیفی با اطلاعات مکانی ادعام گردیده و پیکسل‌ها پراساس شکل. بافت و 
تن خاکستری در سطح تصوير با مقیاس مشخص, قطعه‌بندی شده و طبقه‌بندی تصوير براساس اين قطعه‌ها انجام 
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می‌شود. در قطعهبندى پیکسل‌ها توسط الگوریتم‌های مختلفء با نسبت‌های مختلف طيفى و شكلى طبقه‌بندی و 
براساس ویژگی‌های ab‏ و فضایی در قالب اشيا گوناگون دسته‌بندی شده‌اند. طی اين فرآیند شی‌های تصویری 
متناسب با همگن بودن يا ناهمگن بودن براساس پارامترهای مقیاس» رنگ» شکل» ضریب نرمی و شکل فشردگی 
cal‏ هد at)‏ واد و لا ۱۳۹۶ با آرمون y‏ عطا بهترین انداده که تان هته steal Lob‏ مخقلت go‏ باد 
تعیین شد. در اين پژوهش تصوير قطعه‌بندی شده دارای مقیاس LO‏ شکل ۰/۵ ضریب فشردگی ۰/۵ انجام گردید. 
بعد از مشخص كردن تعداد و کلاس‌های موردنظر در طبقه‌بندی تعيين شی‌های نمونه‌های آموزشی انجام و در 
نرم‌افزار بر سطح تصوير اعمال گردید. سپس طبقه‌بندی با استفاده از نزدیک‌ترین همسایه و روش درخت 
تصميم كيرى صورت گرفت. لازمه استفاده از هر نوع اطلاعات موضوعی, آگاهی از صحت و درستی OT‏ است. دقت 
oS Gre dls sees‏ هل و شیب jd EA‏ 

فرآیند تحلیل شبکه‌ای (ANP)‏ 

روش (ANP)‏ شکل توسعه یافته‌ای از روش (AHP)‏ است که قادر است همبستگی و بازخوردهای موجود بين 
pole‏ در يك تصمیم‌گیری را مدل‌سازی نموده و تمامی تأثیرات درونی اجزای موثر در تصمیم‌گیری را منظور و وارد 
محاسبات نماید. در مدل AHP‏ تنها معیارها با معیار بالایی خود ارتباط داشتند؛ یعنی معیارها به هدفهاء زیرمعیارها 
به معیارها و گزینه‌ها هم به زیرمعیارها وابسته بودند» حتی خود معیارها هم به همدیگر وابسته نبودند. ولی در مدل 
ANP‏ نهتنها خوشه‌ها بر عناصر و عناصر بر گزینه. گزینه بر عناصر تأثیر می‌گذارند. Sh‏ حتی pole‏ بر خودشان و 
بر So‏ حوشه‌ها نیز اثرگذارند. به‌طو ركلى» روش (ANP)‏ دارای مراحلی نیز هست که شامل. ۱- تعيين معیارها و 
شاخص‌ها ۲- تعیین روابط و ارتباطات بين عناصر و خوشه‌ها: در اين گام» es‏ تصميمكيرى به ساختار شبکه‌ای 
تجزيه می‌گردد.۳- مقايسات زوجى بين عناصر و خوشه‌ها: ماتریس‌های مقايسات زوجى مشابه با روش (AHP)‏ 
است و زوجهاى عناصر تصميمكيرى در هر خوشه به نسبت اهميتشان در جهت شرطهاى JES‏ آن‌ها مقايسه 
مىشوند. 6- تشكيل سويرماتريسها: الف- تشكيل سويرماتريس اوليه يا غير وزنى ب- تشكيل سويرماتريس وزنى 
ج- تشكيل سويرماتريس حدى: سويرماتريس وزنی» به توان حدی مىرسد تا عناصر ماتريس همگرا شده و مقادير 
سطرى آن با هم برابر شوند. ۵- نتايج ماتريس خوشهها. T‏ نتیجۂ نهایی عناصر: ضرايب سويرماتريس در ضرايب 
ماتریس خوشهها نرمال شده و درنهایت. نتيجة نهايى عناصر و اولويت آن‌ها مشخص می‌گردد (فرجى سبكبار و 
همکاران ۱۳۸۹). 

گام‌های روش ماباک (MABAC)‏ 


۱- تعيين معیارها و گزینه‌های پژوهش. Y‏ تشکیل ماتریس تصمیم: ماتریس تصمیم در اين روش به‌صورت 
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سلول نيز در واقع امتياز هر كزينه نسبت به هر معيار است. ۳- در اين مرحله نسبت به استاندارسازی و يا 
به‌عبارت‌دیگر بىمقياس نمودن معيارهاء اقدام م ىكردد. در اين يزوهشء از روش فازى و تابع عضويت خطى 
(Linear)‏ و سيكموئيدال (Sigmoidal)‏ به‌منظور استانداردسازى استفاده شده است. در مجموعه‌های فازی, مقدار 
يك به‌عنوان بیشترین ارزش به حداکثر عضویت و صفر به‌عنوان کمترین ارزش» به حداقل عضویت لحاظ می‌گردد 
(حسینی و همکاران (MA‏ 

6- وزن‌دار كردن ماتریس نرمال: در اين گام با استفاده از رابطه زیر ماتریس نرمال را وزن‌دار مىكنيم (رابطه AY‏ 
در مطالعه حاضر جهت وزن‌دهی معیارها از روش ANP‏ استفاده شده است و نتایج حاصل از وزن‌دهی ANP‏ با 
توجه به رابطه زین به تک تک لایه‌های اطلاعاتی لحاظ گردیده است. 

(V_ij=W_j*(nij + 1 ) ١ رابطه‎ 

كه در آن nij‏ عناصر ماتريس نرمال شده N)‏ را Wiss p OUS‏ ضرايب وزنی معيار را OUS‏ می‌دهد. 

۵- تعيين مرز ناحيه شباهت ماتريس (6): در رابطه vijé‏ عناصر ماتريس وزنى (V)‏ و m‏ تعداد كل گزینه‌های 
جايكزين را نشان می‌دهد. پس از محاسبه مقدار gi‏ ماتريس محدودههاى تقريبى مرزى با توجه به معيارهاى LG‏ 
فرمت ۸*1 شكل م ىكيرد و 7 تعداد كل معيارهايى را كه گزینه‌های ارائه شده برای آن‌ها انتخاب می‌شود را OUS‏ 
می‌دهد. 

رابطه G-D^mv ij (Vm) Y‏ ]120 ع) 

بنابراين اگر 2 معيار داشته باشيد يك ماتريس Glxn‏ به‌صورت زیر خواهيم داشت (رابطه ۳): 


[G-Igl.g2.....gn Y رابطه‎ 

aal uale‏ گر ها قاس اھک a‏ ادم شش با o las ala: ada Sana‏ اسرد ور 
دست می‌آید. در واقع بايد ماتريس وزندار را از ماتريس 8 كم کرد. 

Q-V-G ٤ رابطه‎ 

(GI وه بای ات‎ Laos, sss Ale sl انا‎ Oo e Dae sess Jos 
وضعيت هر گزینه را مشخص كرد.‎ 

۷- انتخاب گزینه بهینه: در روش ماباک MABAC‏ با استفاده از رابطه ف امتیاز نهایی هر گزینه را مشخص کرده 
تس کرت وای ی ر 

S_FZ(qij);=1.2....n0j=1.2....m ۵ abl, 


سال يازدهم ارزيابى و يهنهبندى خطر فرونشست با استفاده از .... oy‏ 


محاسبه مقادير توابع kaa‏ توسط كزينهها به‌عنوان مجموع فاصله‌های جايكزين از مناطق تقريبى مرزى 41 به 
رابطه n‏ تعداد معیارها را نشان می‌دهد و m‏ تعداد گزینه است (عليزاده و همكاران (Sayi‏ 


جدول ۲- نوع توابع فازی‌سازی لایه‌های موضوعی تأثیرگذار بر خطر فرونشست دشت اردبیل 


متغیر توضیحات نوع alt‏ فازی 
كيت کلاسه‌بندی مجدد و اختصاص کدها بر اساس اهمیت آن در وقوع خطر فرونشست و سپس فازی کردن. | تابع‌بزرگ فازی 
کاربری اراضى اختصاص كدها بر اساس اهميت آن در وقوع خطر فرونشست و سپس فازى كردن. تابع‌بزرگ فازى 
لیتولوژی اختصاص كدها بر اساس اهميت آن در وقوع خطر فرونشست و سپس فازى كردن. تابع Soy‏ فازى 
فاصله از asl yl‏ در مناطق مجاور رودخانه. پتانسیل وقوع مخاطره فرونشست کاهش می‌یابد. تابع نزدیک فازی 
فاصله از گسل در مناطق مجاور گسل. پتانسیل وقوع مخاطره فرونشست افزایش می‌یابد. تابع نزدیک فازی 
فاصله از شهر و 
5 در مناطق مجاور شهر و روستاء پتانسیل وقوع مخاطره فرونشست افزايش می‌یابد. تابع نزدیک فازی 
زوس 
افت سطح آب با كاهش ميزان سطح آب جاههاء خطر فرونشست افزايش مى يابد. ati‏ خطى کاهشی 


معيارهاى مطرح در يهنهبندى فرونشست دشت اردبيل 

ليتولوزى: قديمىترين سازندهاى شناخته شده در دشت اردبيل دولوميت و دولوميت Sal‏ متعلق به ژوراسیک 
است. رخساره‌های کرتاسه» ماسه سنگ و كنكلومراى قرمز می‌باشد. رخساره‌های دوران سوم شامل ماسه‌سنگ‌های 
توفی و توفی برشى با ميان لایه‌های گدازه آندزیتی» آهک كريستاليزهه آندزیت بازالت با پیروکسن زياد مربوط به 
پالئوژن می‌باشد. آبرفت‌های کواترنر که شامل تراس‌های قدیمی» رسوبات جوان رودخانه‌ای» رسوبات آهکی تراورتن» 
توف» ASE‏ لاهار» مارن و رسوبات جوان دشت بیشترین سهم را در دشت اردبیل و بخش مرکزی أن که بیش‌ترین 
ميزان تمرکز چاه‌ها و قنات‌ها را دارده از رسوبات آبرفتی تشکیل شده است كه در مرکز و قسمت‌شرقی دشت 
بيشترين میزان فرونشست دیده می‌شود و در ساير بخش‌ها به‌واسطه شیب بیش‌تر و نبود رسوبات آبرفتی زمینه برای 
وقوع فرونشست بسیار کمتر از ساير بخش‌ها است (شکل ۲ الف». 

شیب: دشت اردبیل در شمال‌غرب ايران و در قسمت‌شرقی فلات آذربایجان به‌صورت حوضه بسته‌ای است که از 
شمال و غرب منطقه به كوه سبلان و کوه‌های اطراف أن و از شرق و جنوب به رشته‌کوه‌های تالش منتهی می‌گردد؛ 
بنابراین زمین‌های حاشیه دشت از اطراف به سمت دشت شیب دارند. بررسی وضعیت ارتفاع OUS‏ می‌دهد که در 


بخش مرکزی دشت کم‌ترین ميزان ارتفاع وجود دارد. شواهد فرو نيز بیش‌تر در اين طبقه ارتفاعی دیده 


می‌شود. در مورد لايه شیب. محدوده شیب ۰-۵ کمترین میزان شيب در محدوده منطقه می‌باشد که اين شيب نقش 


oí‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


مؤثرى در جمع‌آوری منابع آب زيرزمينى فراهم می‌آورد. هر جه ميزان شيب افزايش می‌یابد. از ميزان تمركز چاه‌ها 
كاسته و ميزان فرونشست زمين هم کاهش می‌یابد (شکل Y‏ ب). 

فاصله از رودخانه: عامل نزديك به رودخانه به‌عنوان عامل مزاحم در ايجاد فرونشست شناسایی كرديد (شايان و 
همکاران. ۱۳۹۵) و با فاصله گرفتن از رود امتياز بیشتری به زمينهاى منطقه تعلق م ىكيرد. هر جه از آبراهه‌ها فاصله 
مى كيريم نياز استفاده از آب‌های زيرزمينى به دليل عدم وجود و استفاده از آب‌های سطحی بالا مىرود و در بى 
استفاده نادرست از آب‌های زيرزمينى وافت سطح آبهاى زیرزمینی» ميزان احتمال فرونشست ففزايش مى يابد 
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شکل ۲- نقشه معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل؛ الف) نقشه لیتولوژی. ب) نقشه شیب. ج) 


نقشه آبراهه 


فاصله از شهر و روستا: عوامل انسان‌ساحت مانند زمین‌های کشاورزی ابنیه بزرگ و صنعتی» تأسیسات زیرزمینی 
مانند تونل» خطوط انتقال نیرو و عملیات حفاری و معدن کاری» سکونت‌گاه‌های شهری و روستایی و شبکه‌های 
حمل‌ونقل جاده‌ای با افزايش مصرف آب‌های سطحیی زیرزمینی» تحمیل بار وارده بر زمين در اثر احداث 
ساختمان‌های بزرگ و تراکم زیرساخت‌ها میزان آسیب‌پذیری اين واحد را افزايش داده و سبب تشدید فرونشست 
شده‌اند (شکل ۳ الف). 

بررسی تراکم چاه‌ها و افت سطح آب: براساس نتایج ارائه شده از مطالعات محققان, به نظر می‌رسد که اعمال 


فشار بي شازحد بر آبخوان ممکن است موجب بروز پدیدۀ فرونشست شود. ازاین‌رو به بررسی وضعیت سفره آب 


سال يازدهم ارزيابى و يهنهبندى خطر فرونشست با استفاده از .... 00 


زيرزمينى در محدودهٌ موردمطالعه پرداخته شد. بهمنظور بررسى و ارزيابى نوسانات سطح آب زیرزمینی» از آمار 
چاه‌های بيزومترى موجود در دشست. ميزان افت سطح آب زيرزمينى طى سالهاى WV‏ تا ۱۳۹۹ محاسبه شد. طبق 
بررسی‌های انجام شده مشخص است که سطح col‏ دشت. رو به پایین رفتن است. به‌طوری‌که کمترین افت سطح 
آب ۱/۳ متر و بيشترين ۲۳ متر می‌باشد. علت اين افت» برداشت‌های غیراصولی به دلیل افزايش حفر چاه‌های عمیق 
و نیمه‌عمیق به‌صورت مجاز و غیرمجاز: کاهش بارندگی در سال‌های اخير و به‌تبع آن کاهش آب‌های سطحی در 
محدوده موردمطالعه می‌باشد. براین‌اساس بیش‌ترین افت در بخش‌های مرکزی و شرقی دشت مشاهده می‌شود که 
احتمال وقوع فرونشست در این قسمت از دشت افزايش مىيابد (شکل ۲ ب). 

فاصله از گسل: به دلیل اثر گذاشتن فعالیت گسل‌ها (نظیر زلزله» بر روی فرونشست زمین» بررسی گسل‌های 
منطقه ضرورت می‌یابد. گسل‌ها به نوبه خود نقش مؤثر و مهمی در ميزان توسعه فرونشست در سطح دشت دارد. 
به‌گونه‌ای که هر جه میزان فاصله از حطوط گسل بیش تر باشد نشان‌دهنده اين است که منطقه تناسب کمتری برای 
احتمال فرونشست زمين را دارا بوده و هر جه میزان فاصله به خطوط گسل نزدیک‌تر شود احتمال وقوع فرونشست 
بیش‌تر می‌شود. در دشت اردبیل نيز به علت وجود گسل‌های فراوانی نظیر اردبیل» بزغوش, آستارا و.؛ که دارای روند 


شمال خاورى - جنوب باختری می‌باشد و به همین دلیل ضریب آسیب‌پذیری بالایی دارد (شکل ۳ ج). 
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شکل -Y‏ نقشه معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل؛ الف) نقشه شهر و روستا ب) نقشه افت 


o‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


کاربری اراضى: انواع كاربرى در دشت اردبيل به‌صورت نواحى انسان‌ساخت» مرتع خوب. مرتع متراكم» مرتع 
متوسطء مرتع ضعیف» پوشش برف دائمی. دیم‌زار. مناطق کشاورزی» زمين باير» باغات و پوشش جنگلی» شهرک 
صنعتی» مناطق آبى و اراضى pl‏ می‌باشد. ازآنجهت دشت اردبيل داراى توي وكرافى هموار می‌باشد و از آبرفت‌های 
جوان و تراسهاى قديم و جديد در اثر ته‌نشست آبرفت‌های رودخانه‌ای با شيب كم تشكيل شده است محل مناسبى 
براى فعالیت‌های زراعى و باغات ايجاد كرده است. لذا اغلب اراضى كشاورزى و نيز باغات در اين واحد گسترده 
شده است. بررسى معيار كاربرى اراضى SE‏ از اين موضوع است كه کاربری اراضى و مراتع و نواحى انسان‌ساز و 
شهرک صنعتی بيش ترين درصد از مساحت مناطق پرخطر را به خود اختصاص داده است؛ كه با مصرف بالاى منابع 


آب زيرزمينى جهت مصارف روزانه و آبیاری احتمال وقوع فرونشست را افزايش مىيابد. 
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شکل ۶- نقشه معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبيل؛ نقشه كاربرى اراضى 


۳- نتایج و بحث 


ازاین‌ری لازم است معیارها قبل از ترکیب با یکدیگر استاندارد شود. در اين مطالعه. مرحله مربوط به ارزش گذاری و 


استانداردسازی به‌صورت توأم و بر مبنای ارزش عضویت در مجموعه فازی در نظر گرفته شده است. در پژوهش 


سال يازدهم ارزيابى و يهنهبندى خطر فرونشست با استفاده از .... ov‏ 


حاضس ارزش عضويت بين مقادير ١ Ge‏ قرار دارد. نقشه‌های استاندارد شده مربوط به معيارهاى مطرح در مطالعه 
حاضر در شكل ۵ نشان داده شده است. 
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aT‏ لاريريى اراضى نفشه اسدابدارد شده اسل از شهر و روسنا 
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نقشه استاتدارد شده فاصله از سل 


LX AA 
LIS 


w^ 


minm 


wu 
wie 


mows avers Amir we ee eee moms | 


شكل ۵- نقشه‌های استاندارد شدۀ فازى معيارهاى مطرح در پهنه‌بندی فرونشست دشت اردبيل 


0۸ جغرافیا و مخاطرات محیطی ob‏ چهارم 


غيروزنى يا اوليه مربوط به مدل» سويرماتريس غيروزنى و سويرماتريس حدی» آخرين مرحله برای تعيين ارزش و 
ضريب نهايى عناصرء محاسبة نتايج ماتريس خوشه‌ها و نرمال‌سازی ضريب عناصر در سويرماتريس حدى توسط 
ضریب خوشه‌ها می‌باشد. نتیجۀ مقایسات زوجی بين گروه‌ها که حدول اوزان گروه‌ها نامیده می‌شود. در جدول )« 


آورده نله Cael‏ 


جدول ۳- ماتریس اوزان گروه‌ها 


خوشه‌ها انسانی زمین‌ساختی فیزیوگرافی هیدرولوژی 

۳۲ vvv Zn ۰,۲ انسانی‎ 
oe D OM SPA زمين ساختى‎ 
تليق‎ efe AN JANE AS فيزي و‎ 
۳۲ ۳۳۳ Ur 0 هيدرولوزى‎ 


نتيجة نهایی عناصر: ضرايب سويرماتريس در ضرايب ماتريس خوشهها نرمال شده و درنهایت. نتيجة نهایی 
عناصر و اولويت آن‌ها مشخص شد (جدول (E‏ كه معيارهاى افت سطح آب چاه‌های موجود در منطقه» فاصله از 


رودخانه» ليتولوزى و كاربرى اراضى به ترتيب بيشترين تأثير را در فرونشست دشت اردبيل را داشتهاند. 


جدول 5- وزن نهايى حاصل از روش ANP‏ 


معیار فاصله از شهر و روستا | کاربری اراضی فاصله از گسل لیتولوژی شيب | افت سطح آب | فاصله از رودخانه 


وزن نهايى ۷ 4۷ (Y‏ 10۲ ۰/۱ ۰/۲۳۱ 2۳۹ 


يس از اعمال كردن وزن‌های حاصل از روش ANP‏ در هر يك از نقشه‌های استاندارد شده و پیاده‌سازی ديكر 
مراحل مدل MABAC‏ نقشه نهايى بهنهبندى فرونشست در ۵ Ab‏ بسيار يرخطر تا بسيار كم خطرء به دست آمد 
(شکل ۱). با توجه به نتايج به دست آمده 6 كيلومترمربع از مساحت شهرستان در Abb‏ بسيار پرخطر و 
۵۹ کیلومترمربع. در طبقه پرخطر قرار دارد. همچنین ۳۸۳/۷۵ کیلومترمربع از مساحت دشت. دارای خطر متوسط 
و به ترتیب 14۱/۲۳ و ۲۸۱/۷۱ کیلومترمربع دارای خطر بسیار کم و کم می‌باشند (جدول ۵). 
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شکل -١‏ نقشه پهنه‌بندی Ja‏ فرونشست دشت اردبیل 


جدول ۵- مساحت مربوط به هر یک از طبقات خطر 


طبقه خطر بسیار پرخطر پرخطر خطرمتوسط كم خطر بسیا ركم خطر 
مساحت به کیلومترمریع ۱۹ ۳۷/۰/۵۹ ۳۸۵ AW TIS‏ 
مسالحت به کر vv Yo/ov ۲۳۳/۸ 1/4۷ VEA‏ 


طبق بررسىهاى انجام شده و نتايج به دست آمده ۳۵/۵۵ درصد از منطقه موردمطالعه در يهنه بسيار يرخطر و 
يرخطر قرار دارد. با بررسى و محاسبه موردى معيارهاى تأثي ركذار بر فرونشست و يراكنش آن‌ها در يهنه بسيار يرخطر 
می‌توان بیان كرد كه فاصله از رودخانه كه طبق جدول (E)‏ مهمترين و تأثیرگذارترین معيار در فرونشست منطقه 
موردمطالعه مىباشد و هر جه از آبراهه‌ها فاصله می‌گيريم نياز استفاده از آب‌های زيرزمينى به دليل عدم‌وجود و 
استفاده از آب‌های سطحى بالا مىرود و در بى استفاده نادرست از آب‌های زيرزمينى و افت سطح آب‌هاء ميزان 
احتمال فرونشست افزایش می‌بابد. لذا ۶۲/۳۲ کیلومتر از منطقه بسیار پرخطر در فاصله ۱۸۰۰۰ متری از رودخانه, 
۸ کبومتر در فاصله ۰ متری و ۱۱/۷۱ کیلومترمربع در فاصله ۰ مترى از رودخانه قرار دارد. هر جه 
میزان فاصله به حطوط كسلى نزدیک‌تر باشد احتمال وقوع فرونشست افزايش می‌یابد و هر جه میزان فاصله از 
خطوط گسلی بیشتر باشد نشان‌دهنده اين است که منطقه تناسب کم‌تری را برای احتمال فرونشست دارا می‌باشد. 
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۱ كيلومتر از يهنه بسيار يرخطر در فاصله ۵۰۰۰ مترى از گسل» ۸۰/۹۸ كيلومتر در فاصله ۱۰۰۰۰ متری»۲۱/۳۵ 
كيلومتر در فاصله ۱۵۰۰۰ مترى و ٠/86‏ كيلومتر در فاصله ۲۰۰۰۰ مترى از كسل قرار دارد. بيشترين احتمال وقوع 
فرونشست در محدوده شیب بين ٠‏ تا ۵ درصد اتفاق مىافتد زيرا اين شيب نقش مؤثرى در جمع‌آوری منابع آب 
زیرزمینی فراهم می‌آورد. هر جه میزان شیب افزایش می‌یابد از میزان تمرکز چاه‌ها کاسته و میزان فرونشست زمین هم 
کاهش مىيابد. ۱۲۷ کیلومتر از منطقه بسیار پرخطر در محدوده شیب ٠‏ تا ۱۵ درصدء ۸۱/۹ کیلومتر در شیب ۱۵ تا 
۰ درصد و ۳۵/۱ کیلومتر در محدوده شیب ۳۰ تا £0 درصد قرار دارد. عوامل انسان‌ساخت مانند ابنیه‌های بزرگ و 
صنعتی. تأسیسات زیرزمینی» عملیات حفاری, شبکه‌های حمل ونقل جاده‌ای و سکونت‌گاه‌های شهری و روستایی با 
تحمیل بار وارده بر زمین میزان آسیب‌پذیری را افزايش داده و سبب تشدید فرونشسست می‌شود به همین دليل 
بيش ترين میزان فرونشست در پهنه بسیار پرخطر در فاصله ۲۰۰۰۰ متری از شهر و روستا با ۱۳۶/۰۹ کیلومتر در 
فاصله 2۰۰۰۰ متری از ۸۰/۲۰ کیلومتن در فاصله ۰۰۰۰ متری ۲۱ کیلومتر و در فاصله ۸۰۰۰۰ متری 51 کیلومتر 
قرار دارد. آبرفت‌های کواترنر كه شامل تراس‌های قدیمی. رسوبات جوان رودخانه‌ای» رسوبات Sal‏ تراورتن» توف 
خاكسترء لاهار, مارن و رسوبات جوان دشت بیشترین سهم را در دشت اردبیل و بخش مرکزی آن که بیش‌ترین میزان 
تمرکز چاه‌ها و قنات‌ها را دارده از رسوبات آبرفتی تشکیل شده است که در مرکز و قسمت‌شرقی دشت بیش‌ترین 


میزان فرونشست دیده می‌شود و در سایر بخش‌ها به‌واسطه شیب بیش‌تر و نبود رسوبات آبرفتی زمینه برای وقوع 


TC‏ 9 بسیار کم‌تر از ساير بخش‌ها است. رسوبات جوان رودخانه‌ای و دشت با ۲۰ کیلومتر. مارن. 
توفهاء سنگآهک و مارن و ماسه‌سنگ توفى ۳۲ كيلومتر از مساحت پهنه بسيار يرخطر را شامل می‌شوند. 

بررسی معیار کاربری اراضی SE‏ از اين موضوع است که مناطق کشاورزی WAL‏ کیلومتر نواحی انسان‌ساز 
۹ کیلومتر دیم‌زار MIT‏ کیلومتر مناطق آبی ۹۸/۸ كيلومترء مراتع خوب ۱۵۵/0 کیلومتر» مراتع متوسط ۱۷۸/۸ 
شهرک صنعتی ۱۳۲/۹ کیلومتر مساحت مناطق بسیار پرخطر تا کم daz‏ را به خود اختصاص داده است. مناطق دارای 
کاربری اراضی کشاورزی آبی به دلیل بهره‌برداری از مناطق آب زیرزمینی» پتانسیل بالایی جهت وقوع فرونشست 
دارد. همچنین مناطق سکونت‌گاهی نيز به دليل ساخت‌وسازها و وجود ابنیه‌های بزرگ و سنگین مستعد وقوع 
فرونشست هستند. با توجه به ايتكه عمده محصولات کشاورزی زیر کشت دشت اردبیل گندم جو و سیب‌زمینی 
می‌باشد که نياز آبی بالایی برای رشد دارند و کاهش چشم‌گیر بارندگی و کمبود منابع آب‌های سطحی در سال‌های 
اخیر لزوم استفاده از منابع آب‌های زیرزمین را بالا برده و همین مساله سبب افت بیش‌تر منابع آب زیرزمینی در منطقه 


شده است. نتایج کاربری اراضی محدوده موردمطالعه و پراکنش فرونشست در اين پژوهش OUS‏ می‌دهد که به دلیل 
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رشد سريع جمعيت و توسعه صنعت و كشاورزى» تقاضای آب نیز بسيار زياد می‌باشد. بااین‌حال بهرهبردارى 
بي شازحد آب‌های زيرزمينى به دليل كمبود منابع آب‌های سطحى آشکار است كه منجر به فرونشست زمين می‌شود 
و باعث آسيب به جادههاء يلهاء لوله‌های آب زيرزمينى و ساير تأسيسات شهرى می‌شود. 


شكل 7- تغييرات ميانكين سطح ايستابى چاه‌های پیزومتری دشت اردبيل 


در دشت اردبیل قبل از دهه ٠١‏ و ۷۰ تعداد چاه‌های عميق و نيمهعميق در مقايسه با چاه‌های دستى يا نیمه عميق 
انگشت‌شمار بود. درنتيجه برداشت از منابع آب زيرزمينى در حد محدود انجام می‌گرفته است. از سال 7/8 به بعد 
استفاده از وسایل و دستگاه‌های حفاری, تعداد چاه‌های حفر شده در دشت اردبیل رو به افزايش نهاد و با رونق 
کشاورزی از دهه ۷۰ به بعد و به‌تبع آن افزایش بی‌رویه برداشت از سفره مذکور شروع به افت سطح آب زیرزمینی 
کرد. تداوم اين وضعیت در سال‌های بعدی به‌حصوص از سال ۱۳۸۵ به بعد باعث بحرانی شدن وضعیت اين دشت 
از بابت افت سطح آب و فرونشست شده است. کاهش آبدهی caola‏ چشمه‌ها و قنات‌ها و تغيير محل و کف‌شکنی 
چاه‌ها و درنتیجه تحمیل هزینه‌ای گزاف در اين خصوص به صاحبان چاه‌هاء افت کیفی منابع آب و شور شدن 
تدریجیی ایجاد محدودیت در مصارف کشاورزی» صنعت و شرب نشست‌زمین و کاهش ظرفیت حجم مخزن 
طبیعی از تبعات افت سطح آب‌های زیرزمینی در دشت اردبیل است. با استفاده از اطلاعات به دست آمده می‌توان 
كفت که قسمت شرقی و مرکزی دشت اردبیل بيش ترين تغیبرات را به لحاظ افت سطح آب زیرزمینی داشته‌اند كه به 
علت مساعد بودن شرايط کشاورزی و بالا بودن سطح زیر کشت محصولات با نیاز آبی بالا و تراکم چاه‌های عمیق 
حفر شده برداشت از آب زیرزمینی در منطقه ذکر شده فراوان بوده و به همین دلیل میزان افت سطح آب‌های 


زیرزمینی بالا می‌باشد. جهت بررسی وضعیت و تأثیر افت سطح آب‌های زیرزمینی و تأثیر آن بر میزان فرونشست 
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دشت اردبيل از اطلاعات ۳۹ چاه پیزومتری (اخذ شده از سازمان آب منطقهاى اردبیل)» طى سال‌های 51 تا 94 (شكل 
۷ كه نتايج تحليل از كاهش £A‏ درصدی سطح ایستابی را در سال ٩‏ نسبت به سال QUA TM‏ می‌دهد. بيشينه سطح 
ایستابی چاه‌های پیزومتری از ٤۵‏ متر به بیش از ۷۰ متر در سال ۹٩‏ رسیده است که حاکی از وخامت اوضاع آبخوان 
دشت اردبیل است. 

بررسی‌های زمینی مؤيد وجود تعداد کثیری شواهد و رخنمون‌های مورفولوژیکی ناشی از افت سطح آب‌های 
زیرزمینی و نشست زمین در دشت اردبیل می‌باشد. ازجمله می‌توان به درز و ترک‌های روی دیوار منازل (شکل ۸) و 
ترک و نشست جاده‌ها و خیابان‌های شهر اردبیل اشاره کرد و ترک‌هایی که در روی زمین‌های کشاورزی شکل گرفته 
سبب توسعه فرسایش شده است. به‌گونه‌ای که اين درز و ترک‌ها باعث کاهش کیفیت زمین‌های کشاورزی و از بين 


رفتن آن‌ها و آسیب به سیستم آبیاری شده است. 


که“ 


شكل ۸- ترک‌های شكل كرفته ناشی از فرونشست در دشت اردبيل 


ci‏ جمع بندى 

در اين پژوهش خطر فرونشست در دشت اردبیل» با استفاده از تلفيق روش‌های ANP‏ و MABAC‏ موردبررسى 
قرار كرفته است. نتايج مطالعه نشان colo‏ عوامل افت سطح آب (جزو عوامل هيدرو اقلیمی)» فاصله از آبراهه (جزو 
خوشه عوامل هيدرو اقلیمی» لیتولوژی (جزو خوشه زمین‌ساختی» به ترتيب با مقدار ۰/۲۲۱ ۰/۱۹۱ ۰/۱۵۲ 
پیشسترین ضريب وزنى را به خود احتصاص دادند. فرزین LS‏ و همکاران (۱۶۰۰) نيز در مطالعه‌ای که در آن به 
پهنه‌بندی خطر فرونشست در دشت جوین با استفاده از مدل تحلیل شبکه‌ای — فازی پرداخته‌اند؛ سطح جاه‌های 
پیزومتری» تراکم چاه‌ها و فاصله از رودخانه را به‌عنوان مهم‌ترین عوامل بیان کرده‌اند. شایان و همکاران (۱۳۹۵)» در 
مطالعه‌ای با عنوان پهنه‌بندی خطر فرونشست در حوضه‌ای آبخیز قره‌چای همدان به اين نتیجه رسیدند که تغییرات 
ناهمسان در داده‌های مورفومتری شبکه زهکشی مهم‌ترین فاکتور موثر در تشدید افت آب‌های زیرزمینی و وقوع 
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فرونش ست در منطقه بوده و عوامل ليتولوزيكى و انسانی بهتنهايى نقش ندارند. با توجه به بررسى نقش و اهميّت هر 
يك از عوامل دخيل در شكل كيرى فرونشست و نيز نتایج حاصل از بهكاركيرى روش تحقیق و تجزيهوتحليل نهايى 
با استفاده از روش ماباک بدعنوان یکی از روش‌های تحليل تصمیم‌گیری چندمعیاره به ترتیب ۲٤٤/۲۹‏ و ۳۷۰/۵۹ 
کیلومترمربع از محدوده مطالعاتی در طبقات بسیار پرخطر و پرخطر قرار دارد. با توجه به نقشة پهنه‌بندی فرونشست. 
دشت اردبیل دارای احتمال Jos‏ باه مناطق بسیار پرخط به‌طور عمده در شیب ۰ تا ۱۵ درصد. سازندهای نامقاوم و 
نیمه‌مقاوم. كه Gase‏ سازندهای سست کواترنری (رسوبات جوان رودخانه‌ای, رسوبات Sal‏ تراورتن. توف اهار 
مارن)» کاربری‌های زراعی: انسان‌ساخت و مراتع» محدوده چاه‌های پیزومتری و Grok‏ محدوده نزدیک گسل و در 
مناطق با فواصل دور از آبراهه‌ها قرار دارند. همچنین می‌توان كفت که نتایج حاصل از اين مطالعه. حاکی از توان GUL‏ 
sli Bits on rcs‏ تشم UU ath‏ وهای برعظر شین ad‏ در این يوسن اراضی مهف که 
نباید کارهای غير اصولی از قبیل برداشت بی‌رویه آب‌های زیرزمینی. حفر چاه‌های عمیق و نیمه‌عمیق» کاشت 
محصولات با نیاز آبی بالا انجام گیرد. همچنین اين امر نیازمند شناسایی عوامل تأثي ركذار, اقدامات کاهش‌دهنده و 
اقدامات تطبیق دهنده ازجمله محدود يا ممنوعیت استخراج آب‌های زیرزمینی در مناطق در معرض فرونشست و 
اجرای مقررات سخت گیرانه جهت استخراج‌کنندگان منابع آب زیرزمینی در مناطق در معرض فرونشست و پایش 
پیوسته شاخص‌های فرونشست است. به‌طورکلی ایجاد چنین شرایط و وضعیتی در رویداد پدیده فرونشست و 
پیامدهای مخاطره‌آمیز و فاجعه‌بار آن برای كشورء قبل از هر چیزی عزم ملی برای پذیرش و مقابله با اين پدیده 
به‌عنوان یک معضل و مخاطره منتهی به يك فاجعه را طلب می‌نماید. موضوعی که درک أن از یک سو مستلزم 
افزايش مطالعات برای شناسایی کامل نواحی مستعد به‌صورت یک طرح ملی است. از طرفی لازم است زیرساخت 
لازم برای آگاهی‌رسانی و همکانی بودن اطلاع از نتایج و معضلات ناشی از اين مسئله برای مردم در جهت مشارکت 
آن‌ها در جهت جلب مشارکت مردم و مسئولان و تعدیل زیرساخت‌های ایجابی فراهم گردد؛ و نهايتاً برای پیشگیری 
و کاهش و تعدیل حجم دامنه فاجعه محتمل از اين مخاطره. دستگاه‌های ذی‌ربط تمامی مقررات و استانداردهای لازم 
ساخت‌وساز را برای نواحی در معرض فرونشست موردبازنگری قرارداده و اساس را بر شرایط نوين تعریف نماید. 
درنهایت می‌توان اذعان داشت با توجه به اين امر كه فرونشست زمین می‌تواند خسارت‌های جبران‌ناپذیر مالی و جانی 
به دنبال داشته باشد و به بسیاری از سازه‌های سطحی و زیرسطحی در مناطق شهری و حومة OI‏ در محدوده 
موردمطالعه سیب برساند. ازاین‌رو بایستی راهکارهای بلندمدت به‌صورت اصلاح روش‌های مدیریت منابع آب و 
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